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Abstract. The probability of the detection of Earth-hke exoplanets may increase in 
the near future after the launch of the space missions using the transit photometry 
as observation method. By using this technique only the semi-major axis of the 
detected planet can be determined, and there will be no information on the upper 
limit of its orbital eccentricity. However, the orbital eccentricity is a very important 
parameter, not only from a dynamical point of view, since it gives also information 
on the climate and the habitability of the Earth-like planets. In this paper a possible 
procedure is suggested for confining the eccentricity of an exoplanet discovered by 
transit photometry if an already known giant planet orbits also in the system. 
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tection 



1. Introduction 

In the last decade, following the discovery of the first extrasolar planet 
around 51 Pegasi (Mayor and Queloz 1995), more than 168 exoplanets 
have been observed. The detection of exoplanets is of great importance, 
since they form planetary systems around their hosting stars, and by 
studying the main properties of these exoplanetary systems, the char- 
acteristics, formation and evolution of the Solar System can be treated 
and understood from a wider perspective. However, the above picture 
is rather ideal than complete yet, since the exoplanets observed until 
now are mainly Jupiter- like gas giants. This is a consequence of the fact 
that by using radial-velocity measurements, which is the most effective 
ground-based observing technique, it is not possible to detect Earth- 
like planets yet. (We note, however, that the detection of a ~ 8Mq) 
planet by Rivera et al. (2005) indicates that radial velocity searches 
might be able to detect Earth-like planets in the future.) 

The most compelling question related to exoplanetary research is the 
detection of Earth-like planets. Beside their importance in testing and 
improving formation theories of planetary systems, a major question is 
their habitability. If an Earth-like planet orbits in the habitable zone 
of a star, there could be chances of appearing (water based) life on 
its surface. On the habitable zone we mean that region around a star, 
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where liquid water can exist on the surface of a planet, for details see 
in (Kasting et al. 1993). 

In order to find Earth-like planets, there are space missions in con- 
struction and planning phase. Such a mission is COROT (sponsored 
by CNES, ESA and other countries) to be launched in 2006, and the 
Kepler Mission (NASA) with a launch in 2008. These missions will 
use the transit photometry as detection technique, which is based on 
measuring the periodic dimming of a star's light intensity caused by an 
unseen transiting planet. Observations performed by these instruments 
will provide the semi-major axis a of the transiting planet calculated 
from Kepler's third law by measuring the period T of the transits, 
and assuming that the mass m* of the hosting star is known. Since 
rri-t, is known only with limited accuracy, an uncertainty appears in a. 
(We note that can be determined by spectroscopic observations 
and by stellar model calculations.) The error appearing in a and in the 
eccentricity e of the transiting planet due to the uncertainties in stellar 
mass will be estimated in 3.2. 

In this paper we present a procedure which can be used to estimate 
the orbital eccentricity and inclination of a transiting planet, if (i) 
we can measure the duration of the transit, and (ii) there is another 
(giant) planet in the system. We shall see that an equation can be 
derived, which connects the mass and the radius R of the star, the 
semi-major axis a, the eccentricity e, the argument of the periastron 
cj, the inclination i of the transiting planet, and the duration r of the 
transit. In this equation there are three unknowns, e, uo, and i. By 
fixing i, the corresponding (w, e) pairs can be visualized as curves on 
the uj — e parameter plane. Thus the problem is undcrdctcrmined and 
there is no way to give an upper limit for e. We note that during one 
revolution of the transiting planet, in principle there are two minima in 
the light curve of the star. The first minimum occurs when the planet 
moves in front of the star, the second appears when the star is between 
the observer and the planet. The maximum of the light curve can be 
observed when both the star and the planet are visible. By using the 
second minimum observation, the eccentricity of the transiting planet 
can be estimated. However, the maximum and the second minimum 
of the light curve can not be seen in the case of an Earth-like planet, 
since its contribution to the whole light flux is not detectable. (The light 
curve is the variation of the light intensity of the star as a function of 
the orbital phase of the transiting planet.) 

On the other hand, as suggested by planetary formation scenarios, 
we expect that beside Earth-like planets Jupiter-like giant planets can 
also be found in the majority of planetary systems. Having discovered 
an Earth-like planet around a star, by using complementary techniques 
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(as observations by Space Interferometry Mission and ground-based 
Doppler spectroscopy) other more massive planets can be identified 
in the system, and their orbital parameters can be determined too. 
(Although the above scenario sounds very optimistic, it is one of the 
scientific goals of Kepler Mission, see 
http : //kepler . nasa . gov/basis/goals . html.) 

The presence of one or more additional giant planets (beside the 
transiting one) results in that, that both ordered and chaotic regions 
appear in the phase space of the system. If the phase trajectory of the 
Earth-like planet is in an ordered region of the phase space, the motion 
of the planet is stable for arbitrary long times. If the initial conditions 
of its orbit are in a chaotic region of the phase space, the motion of the 
planet can become unstable after some time. In this paper we try to 
exclude those orbital parameters of the transiting planet, which result 
in chaotic motion. We shall demonstrate that in some cases it is possible 
to determine an upper limit for the eccentricity and a lower limit for 
the inclination of the transiting planet. We stress that the eccentricity 
is a very important orbital parameter not only from a dynamical point 
of view, but also in studying the habitability and climatic variations of 
the Earth-like planet. 

The paper is organized as follows: first we derive an equation be- 
tween the duration of the transit and some important parameters of the 
star and the transiting body, then we solve this equation numerically. 
After examining the solutions of this equation, we map the stability 
structure of the system assuming the presence of a giant planet. Then 
we determine lower limits for the inclination and an upper bound for 
the eccentricity of the transiting planet depending on the eccentricity 
and the semi-major axis of the known giant planet. 

2. An equation connecting the parameters of a star, a 
transiting planet and the transit 

In this section we shall derive an equation between the orbital param- 
eters of the transiting planet, the star's mass, and the duration of the 
transit from the geometry of the transit. 

Let us suppose that the star's disc is a circle with radius R, and a 
planet is moving in a front of this disc with an average velocity Vti-. If 
the duration of the transit is denoted by r and the lenght of the path 
of the transiting planet is d (see Figure 1), the following approximation 
holds: 

(1) 
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Figure 1. The transit of a planet in front of the stellar disc. The straight sections 
denoted by i?, p, and d/2 form a pythagorean triangle. 



We note that according to Kepler's second law, the velocity of the 
planet is changing during the transit (except in the case of circular 
orbits), however this change is negligible for moderate eccentricities. 
Prom the triangle in Figure 1 

i?^ = P^+(0'. (2) 

From Equation (2) the lenght of the transit's path d can be expressed 

as 



d = - r2 cos2 i, (3) 

where, according to Figure 2, 

p = rcosi, (4) 

i being the inclination (e.g. the angle between the orbital plane and the 
tangent plane to the celestial sphere), and r the distance between the 
center of the star and the planet at the moment of the middle of the 
transit. 

By using the well known formula for r: 

a(l-e^) ..^ 

■r = -7-, , (5) 

1 + e cos V 

(where a is the semi-major axis, e is the eccentricity, and v is the 

true anomaly of the transiting planet), and Equations (1) and (3), the 
average transiting velocity of the planet {vtr) can be written as 




a(l-e2) 



1 + ecosv 



cos^ i , (6) 
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Figure 2. Side-view of the transit, where r is the distance of the planet from the 
star's center and i is the inclination of its orbital plane. 
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Figure 3. The transit as viewed from above. At the middle of the transit vtr is 
nearly equal to y. The coordinate system rj) is the rotation of the coordinate 

system {x, y) by u). 



where v is the true anomaly of the transiting planet at the middle of 

the transit. 

On the other hand, utr can also be approximated on the basis of the 
two-body problem. In the coordinate system (^,7?), in which the axes 
of the orbital ellipse are on the axes ^ and the components of the 
orbital velocity vector are (Murray and Dermott 1999): 

£ = —\ — sinv , (7) 
VP 

?7 = w — (e + cos vj , 

V p 

where p = a(l — e^) is the parameter of the ellipse and /x = k'^{m^:+mp), 
m^: and nip being the stellar and planetary masses respectively, and k 
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is the Gaussian constant of gravity. Let {x,y) denote the coordinate 
system in which the x-axis is the projection of the line of sight (e.g. 
tlic line connecting the center of the star to the observer) to the orbital 
plane of the transiting planet. From Figure 3 it can be seen that the 
system rj) is just a rotation of the system (x, y) by an angle u, which 
is the argument of the periastron of the transiting planet. Thus in the 
coordinate system (x, y) formulae (7) transform as 

X = cos Lu — f] sin Lo , (8) 

y = ^ sin CO + r) cos co . 

Studying Figure 3 one can find that the average velocity Vfr of the 
transiting planet is almost equal to y, which is the vclociy of the planet 
at the middle of the transit. Then by using Equations (7) and (8) we 
find 

- smvsmo; + . / — (e + cost;) cosw. (9) 
P VP 

Studying again Figure 3, it is also true that at the middle of the transit 
v + u = 360° , (10) 
thus the average orbital velocity of the transiting planet is 

VtT = ^l^i^ + ecoso;) . (11) 

Combining Equations (6), (10), and (11) we obtain a relation be- 
tween the orbital parameters of the transiting planet and the duration 
of the transit: 



(1 + ecosu) \ R"^ 



2-\ 



a(l — e 



1 + e cos LO 



cos2z = 0. (12) 



y a(l-e2)^ ^ T 

This equation has the form 

f{a,e,i,(jj,ii,R,T) = , (13) 

where the unknown quantities are the eccentricity e, the inclination i, 
and the argument of the periastron lo. The other quantities, such as 
the semi-major axis a, the mass parameter {^), the radius of the star 
(R), and the duration of the transit (r) are supposed to be known. 
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3. Analysis of Equation (12) 

In this section we present and analyse the solutions of Equation (12) 
by using a model system. We shall see how the orbital eccentricity e 
of the transiting planet depends on its argument of periastron w, or in 
other words, on the orbital position of the planet during the transit. 
We also study the sensitivity of the solutions to the uncertainties of the 
mass of the hosting star. 

3.1. Solution of Equation (12) 

According to the last paragraph of the previous section, the unknown 
quantities in Equation (12) are the inclination i, the argument of peri- 
astron Lu, and the eccentricity e of the transiting planet. Thus by fixing 
i, Equation (12) can be solved numerically, and the (w, e) pairs of the 
solution can be represented as a curve on the lj — e parameter plane. 

In order to study the solutions of Equation (12), we give specific 
values for the parameters in Equation (12). Let us assume that the 
mass of the transiting planet is 1 Earth-mass, and it revolves around 
a 1 Solar-mass star with radius R = 6.96 x 10^ m, in an elliptic orbit 
characterized by a = 1 AU, e = 0.1, i = 89.95°. We consider two 
cases: (i) oj = 30° and (ii) u = 130°. It can be calculated easily that 
in these cases the durations of the transits are r = 0.488029 day and 
T = 0.563743 day, respectively. 




50 100 150 200 250 300 350 

argument of periastron 



Figure 4- Solutions of Equation (12) for different inclinations when t = 0.488029 
day. The original solution, which results in the above t, is marked with a filled circle 
at a; = 30°, e = 0.1, and i = 89.95°. 

By observing transits caused by the above planet, we can measure 
their duration r and period T, from which the semi-major axis a can be 
calculated. In the first case r = 0.488029 day, and for different values 
of i the corresponding u — e curves are plotted in Figure 4 for e < 0.5. 
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Figure 5. Solutions of Equation (12) for different inclinations when t = 0.563743 
day. The original solution, which results in the above t, is marked with a filled circle 
aXu] = 130°, e = 0.1, and i = 89.95°. 

(We note that it is possible to study larger values of e as well, but for 
our dcmonstrational purpose we restrict ourselves for e < 0.5.) We also 
mark the real (w, e) solution by a filled circle on the curve corresponding 
to i = 89.95°, but as we can see, there is no way to restrict efficiently 
the infinite set of solutions. The only restriction is that the solutions 
can not be chosen from the region above the u — e curve corresponding 
to i = 90°. 

In the second case corresponding to r = 0.563743 day, the u — e 
curves arc plotted in Figure 5. It can be seen that if the transit happens 
nearly at the apastron, the set of the solutions of Equation (12) is more 
limited than in the first case. Only those {to — e) pairs satisfy Equation 
(12), which are in the region above the i = 90° curve. The real solution 
is also marked by a filled circle on the curve i = 89.95°. 

Equation (12) has an infinite set of solutions formed by pairs of {to, e) 
values. If only the duration of the transit is known, it is not possible to 
choose which {u, e) pair represents the real parameters of the transiting 
planet. 

3.2. Sensitivity of the solutions of Equation (12) to the 

STELLAR mass' ERROR 

As mentioned in the Introduction, the mass of the hosting star is known 
only with limited accuracy. This uncertainty in the stellar mass affects 
the semi-major axis of the transiting planet and causes an error in the 
eccentricity estimates. The semi-major axis of the transiting planet can 
be calculated from Kepler's third law 




(14) 
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where T is the period of the transits, is the mass of the hosting 
star, and nip is the mass of the transiting planet, respectively (A; is 
the Gaussian gravitational constant). In the case of Earth- like planets 
nip « m^,, so nip can be neglected. Denoting the errors in a and m* 
by Sa and (5m*, respectively, and expanding Equation (14) around a 
and m* up to first order we obtain: 

3a^6a = 7^S'm^=. (15) 

47r^ 

Dividing Equation (15) by and using again Equation (14) (by ne- 
glecting nip) we find for the relative errors 

3^ = ^. (16) 

a m* 

By using Equations (16) and (12) we calculated numerically the rel- 
ative errors in the eccentricity of the transiting planet as a function 
of oj for the specific case when i = 89.95° and r = 0.488029 day. 
Then we determined the maximum values {Se/e)„iax of the relative 
error in eccentricity for different values of (5m* / m* . We found that for 
(5m*/m* = 3% ((5e/e)max = 12%, for (5m*/m* = 6% ((5e/e)max = 25%, 
and for (5m*/m* = 10% ((5e/e)inax; = 38%. 



4. A possible confinement of the eccentricity of the 
transiting planet 

4.1. The case of one additional giant planet 

In this section we shall investigate the case when, beside the newly 
discovered planet, an already known giant planet also orbits around 
the hosting star. The presence of such a planet makes the problem non- 
integrable and both ordered and chaotic regions can be found in the 
phase space of the system. We suppose that the most probable orbital 
solutions of the transiting planet are those, which emanate from the or- 
dered regions of the phase space. Such values of the orbital parameters 
of the transiting planet which would result in chaotic behaviour are 
unlikely, since in long terms the orbit of the planet could be unstable, 
thus these solutions should be avoided. We expect that the presence 
of a second (giant) planet represents a dynamical constraint, which 
reduces the infinit set of solutions of Equation (12), and gives an upper 
limit for the maximum eccentricity of the transiting planet. We shall 
also demonstrate that by stTidying the solution-curves of Equation (12) 
together with the stability structure of the uj — e plane, a lower bound 



zs_sandor.tex; 5/02/2008; 11:35; p. 9 



10 



Zs. Sandor 



for the inclination can also be determined. In what follows we shall 
investigate the stability in the (w — e) plane within the framework of 
the planar general three-body problem. 

In order to map the stability properties of the (w — e) plane we 
used the Relative Lyapunov Indicator (RLI) (Sandor et al. 2000, 2004). 
The initial lo and e values arc chosen from the intervals e € [0, 0.5] 
and uj G [0°,360°] with Ae = 0.025 and Aw = 2°. The initial value 
of the semi-major axis of the transiting planet is always a = 1 AU, 
while its true anomaly is calculated according to Equation (10) as v = 
360° — uj (see also Figure 3). The masses of the planets are 1 Jupiter- 
mass and 1 Earth-mass, respectively, and the mass of the hosting star 
is 1 Solar-mass. 

For each pair of the initial (oj, e) values we assign the RLI of the 
corresponding orbit calculated for 500 periods of the transiting planet. 
If the RLI is small (^ 10~^^ — 10^^^), the corresponding orbit is ordered 
and stable for very long time. If the RLI ~ 10~^^ — 10^^ the orbit is 
weakly chaotic. In practical sense such an orbit could be (Nekhoroshev) 
stable for very long terms as well, however, it can not be stable for 
arbitrary long time. Thus the regions characterized by these RLI values 
can already be the birth places of unstable orbits. Orbits having larger 
RLI ~ 10^^ — 10^^, are strongly chaotic, and they will be unstable 
after certain time. In our stability maps the ordered regions are denoted 
by light, the weakly chaotic regions by grey, and the strongly chaotic 
regions by dark shades. 

In what follows we consider the cases where the parameters of the 
known giant planet having 1 Jupiter mass are the following: ai = 2.0 
AU, ei = 0.2, and 0.3. We fix the angular elements of the giant planet to 
A = = 0° . In Figures 6 and 7 we show the dynamical structure of the 
uj — e parameter plane for the two different values of ei. In these figures 
we also plot the solution curves of Equation (12) by using r = 0.488029 
day. From Figure 6 it can be seen that there are two upper bounds 
for the eccentricity of the transiting planet depending on whether the 
transit occurs near the periastron, or near the apastron. If the transit 
is near the periastron uj < 80°, the upper limit of the eccentricity is 
e < 0.27, since the uj — e curves cross the chaotic region around this 
value. If the transit happens around the apastron uj G [150°, 220°], the 
upper limit of the transiting planet's eccentricity is e < 0.22. In this 
case a lower limit can also be given for the inclination, i > 89.°85. The 
real solution is marked (as a filled circle) on the curve corresponding 
to i = 89.95°. 

If the eccentricity of the giant planet is ei = 0.3, see Figure 7, the 

maximum upper limit of the transiting planet's eccentricity is e < 0.18. 
However, in this case there exists a lower limit e > 0.05 as well. If the 
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Figure 6. The stability map of the u) — e parameter plane, when ai = 2.0 AU and 
ei = 0.2. The w — e curves for different i are also plotted when r = 0.488029 day. 
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Figure 7. The stability map of the oj — e parameter plane, when a\ = 2.0 AU and 
ei = 0.3. The w — e curves for different i are also plotted when r = 0.488029 day. 



transit took place around the periastron the corresponding oj and e 

values would result in weakly chaotic orbits. A lower bound of the 
inclination in this case is i > 89.89°. Among the two values of the giant 
planet's eccentricity, this latter represents a more effective dynamical 
constraint for the orbital parameters of the transiting planet, which are 
a = 1.0 AU, e = 0.1, u = 30°, and i = 89.95°. 
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Figure 8. The stability map of the u) — e parameter plane, when ai = 2.0 AU and 
ei = 0.2. The w — e curves for different i are also plotted when r = 0.563743 day. 



Next we shall study the cases when the transit takes place around 
the apastron u = 130°, resulting in a transit's duration r = 0.563743 
day. We used the same orbital parameters of the giant planet as in 
the previous case, ai = 2.0 AU, ei = 0.2 and 0.3 respectively. Figure 




orgurnent of pcriostron [deg] 

Figure 9. The stability map of the oj — e parameter plane, when oi = 2.0 AU and 
ei = 0.3. The uj — e curves for different i are also plotted when t = 0.563743 day. 

8 shows the situation when e\ = 0.2. The maximum upper limit of 
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the eccentricity of the transiting planet is e < 0.24 corresponding to 
CO ~ 220°. The lowest bound of its inclination is i > 89.89°. 

Finally, when ei = 0.3, as shown in Figure 9, there are no stable 
solutions for the orbital parameters of the transiting planet. 

We have also investigated cases, when the semi-major axis of the 
giant planet was smaller or larger than 2 AU. If ai is smaller, a smaller 
ei is enough to result in an effective dynamical constraint. If ai is 
larger, the eccentricity of the giant planet should be larger as well for 
an efficient dynamical constraint. 

In the above cases wc assumed that the giant planet and the tran- 
siting planet revolved in the same plane. This, based on Solar System 
examples, seems to be a reasonable assumption, however, the orbital 
planes of the planets may differ slightly from each other. Thus we 
performed mimcrical simulations by using a small mutual inclination 
1 = 5° between the orbital planes as well. Comparing the corresponding 
uj — e parameter planes we have not found any significant differences 
between the planar and the spatial cases with small mutual inclination. 

4.2. The case of two additional giant planets 

Among the extrasolar giant planets with long observational baselines, 
there is a high rate of multiple planet systems. We expect that this rate 
will increase further with the accumulation of the observational data. 
Therefore we studied a case, when beside a transiting Earth-like planet 
there are two additional giant planets. 

The inner giant planet's initial orbital parameters are the same as in 
4.1. The outer giant planet, having 1 Mjup, moves in an orbit charac- 
terized by 02 = 4AU, 62 = 0.05 in the same plane as the transiting and 
the inner giant planet. Figure 10 shows the stability properties of the 
(jj — e parameter plane of the transiting planet. In order to study the 
effect of the outer giant planet on the UJ — e parameter plane, we should 
compare Figure 10 either to Figure 8 or to Figure 6. In these figures 
the initial conditions for the Earth-like and the inner giant planet are 
the same. 

The main structure of these figures is very similar to each other, 
however in the case of the additional outer giant planet the strongly 
chaotic region on the u> — e parameter plane is more enhanced. This 
means that the dynamical constraint on the transiting planet could be 
more efficient in the case of two giant planets. Of course the strongness 
of the dynamical constraint depends on many parameters of the outer 
giant planet. To explore the complete parameter space is beyond the 
scope of this paper, it will be the subject of a forthcoming research. 



zs_saiidor.tex; 5/02/2008; 11:35; p. 13 



14 



Zs. Sandor 




100 200 : 

argument of perlastron [deg] 



Figure 10. The stability map of the u) — e parameter plane, when two giant planets 
are in the system. The semi-major axes and the eccentricities of the giants are 
ai = 2.0 AU, ei = 0.2, 02 = 4.0 AU, and 62 = 0.05. The uj - e curves for different i 
are also plotted when t = 0.488029 day. 

5. Conclusions 

The detection of Earth-like cxtrasolar planets by using ground based 
spectroscopic methods is beyond the present capabilities of observa- 
tional astronomy. In the near future space instruments will be launched 
such as COROT and Kepler Mission which are devoted to observe such 
planets by using transit photometry. 

In this paper we addressed the question whether it is possible to de- 
termine the orbital elements of Earth-like planets discovered by transit 
photometry if, apart from the period, the duration of the transit can 
be measured too. We supposed that the mass and the radius of the 
hosting star is known. We derived an equation, which connects the 
stellar and planetary masses, the duration of the transit, the semi-major 
axis, the eccentricity, the argument of periastron and the inclination of 
the transiting planet. By fixing the inclination, this equation contains 
two unknown variables, the argument of periastron u) and eccentricity 
e of the transiting planet. Thus the solutions for different inclinations 
can be represented as curves on the uj — e parameter plane. 

In the last section of the paper wc assumed that beside the transit- 
ing Earth-like planet additional giant planets orbit around the star as 
well. This assumption is quite reasonable if we accept the formation 
theories of planetary systems supporting the simultaneous presence 
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of both rocky, Earth-like and gaseous, Jupitcr-likc planets. Since the 
detection of giant planets is possible by radial velocity measurements, 
we assumed their orbital parameters to be known. (We note that from 
radial velocity measurements only nip sin ip is known, where nip and 
ip are the mass and inclination of the giant planet. However, based 
on Solar System examples, wc can assume that the giant planets or- 
bit nearly in the same plane as the transiting Earth- like planet, so 
its inclination is ip 90°.) By using the framework of the general 
three-body problem, we investigated the influence of the known giant 
planet(s) on the transiting Earth-like planet's cj — e parameter plane. 
We found that on the lo — e parameter plane beside ordered domains 
chaotic regions appeared as well, where in long terms the motion of 
the transiting Earth-like planet may become unstable. Assuming that 
chaotic behaviour for an observed transiting planet is unlikely, wc could 
determine an upper limit for the eccentricity, and a lower limit for the 
inclination. 

In a future work we plan to extend our studies to investigate system- 
atically the stability structure of the {to — e) parameter plane for various 
values of the giant planet's semi-major axis, eccentricity, and inclina- 
tion. Since the mass of the hosting star and several orbital parameters of 
the giant planets are known only with a limited accuracy from the radial 
velocity observations (see for instance Ford et al. 2005), we also plan 
to follow the propagation of these errors through the method presented 
in this paper. In our future investigations we intend to consider cases 
of more massive transiting planets as well. 
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